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ABSTRACT 

The accretion-powered pulsar 4U 1626-67 experienced a new torque reversal 
at the beginning of 2008, after about 18 years of steadily spinning down. We 
present a spectral analysis of this source using two pointed observations per- 
formed by Suzaku in 2006 March and in 2010 September. We confirm with 
Suzaku the presence of a strong emission- line complex centered on IkeV, with 
the strongest line being the hydrogen-like Ne Lya at 1.025(1.5) keV. We were 
able to resolve this complex with up to eight emission lines. A dramatic increase 
of the equivalent width of the Ne Lya 1.021 keV after the 2008 torque reversal oc- 
curred, reaching almost the same value measured by ASCA in 1993. In addition, 
we confirm the general decrease trend of the equivalent widths during the spin- 
down period. We also report on the detection of a cyclotron line feature centered 
at —37 keV. In spite of the fact that a dramatic increase of the X-ray luminosity 
(0.5-100 keV) of a factor of —3.5 occurred between these two observations, no 
significant change in the energy of the cyclotron line feature was observed. How- 
ever, the intensity of the — IkeV line complex increased by an overall factor of 
- 10 . 

Subject headings: accretion, accretion disks — binaries: — pulsars: individual (4U 1626- 
67 ) — stars: neutron — X-rays: stars 
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1. INTRODUCTION 


4U 1626-67 is a Roche-lobe filling binary syste m, where matter from the companion 
flows from the LI point into an accretion disk ([Reynolds et ah 119971) . This low mass X-ray 
binary (LMXB) consists of a 7.66 s_ X -ray puls ar a ccreting from an extremely low mass 
companion (0.04 M© for i = 18°) (Levine et al. 1 98 81). Although orbital motion has never 
been detected in the X -ray data, pulsed optical emission reprocessed on the surface of 
the secondary reve aled (IMiddleditch et al.l Il98ll) the 42 min orbital period, confirmed by 
Chakrabartvl (119981) . Most likely the b inary system contains a hydrogen -depleted secondary 
to reach such a short orbital period (IPaczvnski and Sienkiewiczl 119811) . The faint optical 
count erpart (KZ TrA, V~17.5) has a strong UV excess and high optical pulse fraction 
( McClintock et ah 1977 . 1980 ). 


4U 1626-67 h as revers ed the sign of the torque in only two occasions since its discovery 
by Uhuru (.Giacconi et al. 119721) . Monitoring of the source by the Burst and Transient 
Source Experiment (BATSE) on board the Compton Gamma Ray Observatory ( CGRO ) 
starting in April 1991, found the pulsar spinning down , implying a changed sign in the 
accretion torque (Wilson et al.lll993 ; Bil dsten et ah 1994 ). It was estimated that the reversal 
occurred in mid-1990. More recently, our daily monitoring of 4U 1626-67 with Fermi/ G BM 
starting in 2008 August, found the pulsar again spinning-up rather than spinning-down. 
Swift / BAT observations allowed us to cover the evolution of this second torque reversal. 
The transition took p la ce at around MJD 54500 (2008 Feb 04) and lasted approximately 150 
days (Camero-Arranz et al. 2010*). 


A 48 mHz quas i -periodic oscillation ( Q PO) has been d e tected in the X-ray emission 
(IShinoda et al.l 11990c IKommers et al.l Il998l) . lOrlandini et al.l (119981) inferred a neutron star 
magnetic field in the range of (2. 4-6.3) xlO 12 G. To compute this magnetic field range a 
source distance of 5-13 kpc was assumed (IChakrabartvl 19981 . and references therein). A 
~37 keV absorp tion cyclotron feature was found in the 0.1-200 keV BeppoSAX spectrum 
(lOrl and ini et al.l 19981 ). The X-ray broad-band continuum was fitted with low-energy ab- 
sorption, a blackbody, a power law and a high energy cutoff. 


Angclini et al,! (19951) using the ASCA observatory discovered for the first time a strong 


emission-line complex centered on 1 keV in the X-ray spectrum of 4U 1626-67. The strongest 
line was identified as hydrogen- like Ne Ly a at 1.008 keV. The strength of the neon emission 
compared to the expected iron L complex implied a large neon overabundance. 


We present a spectral analysis of 4U 1626-67 with Suzaku, using two observations from 
2006 March and 2010 September. This allows us to better characterize the recent torque 
reversal occurred in 2008 February. 
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2. OBSERVATIONS and ANALYSIS 



covering the 10-70 keV band and the gadolinium silicate crystal detector (GSO) covering 
the 60-600 keV band. The XIS has four separate detectors, XIS 0-3, with XIS 1 being a 
backside illuminated CCD and the others frontside. XIS 2 was lost due to a micrometeor 
hit in late 2006. 


The first Suzaku 4U 1626-67 observation occurred in 2006 March (MJD 53797), and 
thus before the torque reversal in 2008. During this observation the spacecraft aimpoint 
was optimized for the XIS detectors, and Suzaku was run in a 1/8 window mode with Is 
integration time, and an effective area of 128 x 1024 pixels. The second observation in 2010 
September (MJD 55444) was also placed on XIS aimpoint but run in a 1/4 window mode, 
with 2 s integration time and an effective area of 256 x 1024 pixels. For this observation only 
3 XIS detectors were available, as we discussed in the previous paragrap h. The Suzaku data 
were reduced with tools from the HEASOFT v6.11 package flArnaudl 1996!) and CALDB 
release 20110804. Following the Suzaku ABC guide]/], the first step of our analysis was to 
reprocess the data by running the FTOOL aepipeline for the XIS, the PIN and the GSO. 
The pipeline performs calibration as well as data screening. 

The attitude of the Suzaku spacecraft exhibits variability over the course of the observa- 
tions and therefore the image of the source is not at a fixed position_on_the CCD. This is due 
to a thermal distorsion of the optical bench. As described in Nowak et ah (2011). standard 
processing reduces this variability and improves the PSF image. However, this processing is 
not enough and thus we used the aeattcor2.sl software. Unfortunately, we could not correct 
fully the double image of the source in the XIS images for the second observation. We used 
the ftool pileest.sl to estimate the amount of pileup in the XIS images and disregarded re- 
gions with >8% pileup fraction during the spectral extraction. For the first observation the 
pileup fraction remained below this limit in the entire source extraction region (a circle of 
~ 70” radius centered on the source image). In the second observation pileup fractions of 
up to 15% were reached in the center of the source image so the extraction region consisted 
of a circle of ~ 130” radius minus a central circular exclusion region with a ~ 35” radius. 

The primary tool for extracting data products (spectra, lightcurves) from XIS data is 
xselect , which is part of the general HEAsoft distribution. Events in the XIS detectors were 


1 http: / /heasarc. nasa.gov/docs / suzaku/analysis /abc 


4 


in the 3x3 and 5x5 editing modes. With xselect we applied to the cleaned events filters 
with the good times intervals. Source and background regions are created from the image. 
The filtered events were then used to create source and background spectra for the desired 
extraction regions. Products from the two editing modes were added. The XIS spectra 
were grouped close to the energy resolution of the detectors (see Nowak e t ah 12 01 IIP . We 


created response matrices and effective area hies with the xisrmfgen and xissimarfgen tools, 
respectively. The calibration of the XIS is uncertain around the SiK edge, thus events in 
the 1. 8-2.0 keV range were excluded in the fitting. 


For the PIN, we extracted spectra from the cleaned event hies using hxdpinxbpi. This 
tool produces the dead time corrected PIN source spectrum as well as the PIN background 
spectrum. The non X-ray background is provided by the Suzaku team and can be down- 
loaded from HEASARC (ftp://legacy.gsfc.nasa.gov, |) archiv^E The PIN spectra were 
grouped to have a signal-to-noise ratio > 10 in each energy bin, and we considered spectra 
between 14-100 keV. We point out that due to the changes in instrumental settings the PIN 
response matrices are different depending on the epoch of the observation. 


The GSO spectra were created starting with the clean event hies using the hxdgsoxbpi 
tool. The background was downloaded from HEASARcjE Response hies were taken from 
the CALDB database and the GSO ARF hie from the teanJE The distributed GSO Non 
X-ray Background is grouped, so GSO spectra have to be grouped accordingly, using the 
spectral binning given by the hie gsobgd64bins.dat available at HEASARcjE GSO spectra 
were restricted to the ~60-100 keV range. 


RESULTS 


The 0.5-100 keV broad-band spectrum obtained using XIS (a combination of detec- 
tors 0 and 3), PIN and GSO was htted in XSPEC 12.7.0. We used a canonical model 


which allowed us t o compare our spectra l study with previous works jAngelini et ah 


r 


Pravdo et al.l Il979t lOrlandini et al.l 119981: iKrauss et al.l 120071: I Jain et al. 


2010 


1995 


among oth- 
ers) and update t he long-term X-ray hux history of 4U 1626-67 relative to the hux mea - 
sured by HEAO 1 ( Chakrabartv et al. 1997 : Krauss et al. 2007 : Camero-Arrauz et ah 2010 k 


2 f/suzaku/data/background/pinnxb_ver2 .CLtuned 
3 t/suzaku/data/backgroimd/gsonxb_ver2.0 

4 http : / / www. astro .isas.ac.jp/ suzaku / analysis/hxd / gsoarf / 

5 f/ docs /suzaku/analysis /gsobgd64bins.dat 
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The first model includes low-energy absorption, a blackbody component, a power law, 
a high-energy cutoff at ~20 keV and a Gaussian absorption line at ~37keV (PHABS 
(GAUSS+BBODY+POWLAW) HIGHECUT*GABS)*CONST. The abundance was set to 
wilm (Wilms et al. 20001 ). A broad line near 6.4 keV improved the fit indicating the presence 
of an iron line. The addition of the 37keV cyclotron line feature was also needed after an 
initial inspection of the residuals (see Fig [I])- However, the y 2 was unacceptable using this 
model. The residuals still indicated a strong line feature at IkeV. Following the study by 


Aneelini et ah (119951 ). we added up to 7 narrow Gaussian lines, with their widths fixed to 


the ASCA values (10 eV for all but the 1.46 keV line with 57 eV), and the only requirement 
for inclusion was reducing the y 2 . We also fit a recombination edge yielding a similar y 2 
compared to that obtained using a broad Gaussian line. The edge (and temperature) ob- 
tained for the first and second observation are 1.45±0.05 keV (42 eV) and 1.47±0.03 keV 
(45^5 eV), respectively. A posteriori inspection of the residuals suggested the presence of 
an extra Gaussian component at energy ~1.34keV. Therefore, up to 8 Gaussian components 
were included in the final fit. The summary with the best fit spectral parameters obtained 
has been divided in two tables for clarity. Table [1] shows the best fit broad band spectral 
parameters and Table [5] summarizes the emission lines detected by Suzaku. As a check, we 
fit the broad-band spectrum using different combination of XIS detectors yielding similar 
results. 


A drastic luminosity increase happened between the two observations, changing from 
L(o. 5— 10 keV) ~2.5xl0 34 d| pc to 10xl0 34 d^ pc erg s _1 . Surprisingly, the value of the center of 
energy of the cyclotron line feature does not exhibit a significant change. On the other hand, 
the N h column density dropped while the black body temperature increased by a factor of 
~ 2. The rest of the broad band parameters remained practically the same before and after 
the 2008 reversal. 


Regarding the line emission at low energies, for the first observation only five of the eight 
fitted lines (seven initially observed by ASCA) were significantly detected. For the ones that 
we could not constraint, they were fixed to the ASCA values (accepting only lower flux 
intensity values). Absorbed fluxes for the first model in the 0.5-10 keV, 2-10 keV, 2-20 keV 
and 2-50 keV bands are 2.08(4), 1.78(2), 3.15(3), 4.81(3) xlO~ 10 erg cm“ 2 s _1 and 1.0(3), 
0.706(9), 1.2(1), 1.37(9) xl0~ 9 erg cm _2 s _1 , for the first and second Suzaku observations, 
respectively. 
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DISCUSSION 


In 1990, the torque reversal experienced by _4U 1626-67 was accompanied by a spectral 
transition and a chang e in the luminosity f Yi and Vishniadll999l) . The p ower law photon in- 


dex changed from 1.6 (IPravdo et al.lll979 


Angelini et ah 


19951) to ~0.6 (Owens et ah 1997). 
The black body temperature also decreased from ~0.6keV to ~0.3keV. Our Suzaku obser- 
vations revealed that after the 2008 torque reversal 4U 1626-67 reached almost the same 
luminosity level as in 1977 (see Fig. [2]), but within a period of only ~2 years. Despite the 
fact that the black body temperature transitioned from ~0.23keV back to ~0.6keV in 2010, 
the photon index of ~1 did not show any significant change. The temperature after the 2008 


torque reversal is in ag reement with the value obtained with RXTE (jCamero-Arrauz et ah 
2010l; I Jain et al.l 120101) . For 4U 1626-67 a photon index of ~1 has been obse rved by sev- 


Camero-Arranz et ah 


Orlandini et al. 

1998 

), Chandra and XMM ( 

Krauss et al. 

2007; 

Schulz et al. 

2001) 


In 


4U 1626-67 during the 2008 torque reversal differed from that expected in models, which sug- 
gested significant changes of the accretion flow structure in spin- up/spin- down transitions. 
The spectrum became hardest during the reversal and the value of the hardness ratio before 
and after these events did not differ significantly. This indicated that the recent torque 
reversal could be associated with changes of physical conditions in the inner part of the disk 
or/and in the region of its interaction with the magnetosphere. 


In the present work, following Angelini et ah| (119951) . we have decomposed the strong 
emission-line complex centered on 1 keV found in the Suzaku X-ray spectrum of 4U 1626- 
67. In both Suzaku observations the strongest line is the hydrogen-like Ne Lya line at 
~1.025keV, in agree ment with the ASCA observations. The broad feature at ~1.48 keV 
found by (Angelini et al.l (119951) was close to the expected energy of the H-like Ne recombi- 
nation edge, H-like and He-like Mg emission at 1.36, 1.47 and 1.34 keV, respectively. The 
measured 57 eV broadening of the line at ~1.48keV also with Suzaku, suggests that we may 
be also detecting the recombination continuum of Ne. Furthermore, with Suzaku an eighth 
line centered at ~1.34keV was detected, which most likely is the He-like Mg line. We would 
like to point out that this line energy is so close to the Ne recombination continuum and the 
Mg K lines, and therefore we cannot exclude a different origin for the detected line. 

During the spin-down period, comparing the strength of the first seven lines detected in 
the first Suzaku observation with previous measurements by ASCA, we see that the overall 
decrease is of the o rder of ^2.5 (see T able [2]) . Be fore the 2008 torque reversal, X-ray mis sions 


like Be yyoSAX (Owens et al.lll997l) . Chandra (ISchulz et al.l l200ll; iKrauss et al.l 120071) and 


XMM ( Krauss et al. 2007h were able to resolve only up to 3 or 4 lines of this complex. 


In 
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general, the measurements of the intensity of those lines are in agreement with our results. 
During the new spin-up period after the 2008 torque reversal, the flux of the lines increased 
by a factor of ~5 except the Ne Ly a line at 1.021 keV that increased by a factor of 10. This 
result is expected taking into account such a change in luminosity. Unfortunately, during 
a previous spin-up period there was only one measurement of the ~1 keV complex by the 
Einstein observatory in 1979, and it was not resolved due to the lower data quality. The 
overall luminosity decrease associated w ith the 1990 torque rever sal was a facto r of in the 


0.5-10 keV range (iAngelini et al.l [19951 . and references therein). Krauss et ah (120071 ) found 
a decrease of the equivalent widths of the NeHea 0.914 keV and NeLya 1.021 keV lines 
during the spin-down period. The behavior of the EW for the O Hea 0.568 keV and O Ly a 
0.653 keV lines was more flat during the same period. With Suzaku we found a dramatic 
increase of the EW of the NeLya 1.021 keV line after the 2008 torque reversal. The EW 
for this line reached almost the same value measured by ASCA in 1993. For the rest of the 
lines the increase is not that pronounced. We also confirm the general decrease trend of the 
equivalent widths during the spin-down period. 

In 1993, the strength of the neon emission in the ASCA spectrum of 4U 1626-67 implied 
a large neon overabundance, sug gesting that the com panion star was burning (or had burned) 
helium. A more recent study by Schulz et al. (_2001) with Chandra found that the emission 
lines had a double-peaked profile. They appear to originate in hot (~ 10 6 K) dense material 
just below the X-ray-heated skin of the outer Keplerian accretion disk, or else possibly in 
a disk wind driven from the pulsar’s magnetopause. Based on the inferred local abundance 
ratios, they claimed that the mass donor is probably the 0.02 M© chemically fractionated 
core of a C-O-Ne or O-Ne-Mg white dwarf which has previously crystallized. 

Spectroscopic observations of this iron K-fluorescence line are one of the mo st us eful 


probes to investigate the matter distribution surrounding the X-ray binary system. Kii et al . 


(119861) concluded that the matter accumulated on the Alfven shell seems to be the most likely 
site to reprocess the iron k-fluorescence line in disk-fed X-ray pulsars like 4U 1626-67, rather 
than the accretion disk and/or the atmosphere of the companion star. They obtained an 
upper limit of EW<60eV for this line. With Suzaku we were able to find positive evidence 
for the iron emission line around 6.4 KeV. The iron K-fluorescence line from neutral or 
l owly ionized iron has been observed i n 4U 1626-67 in the past by different obser vatories 
( Pravdo et al. 1979 : Owens et al. 1997 : Camero-Arranz et ali 2010 : Jain et al. 2010 b In the 
present observations we found that the EW changed from ~70 eV to 24 eV, the opposite to 
the expected behavior (that is, no change). These two different values for the iron line EW 
might suggest that changes in the matter accumulated on the Alfven shell occurred after the 
2008 torque reversal. 


Orlandini et ah (1998) discovered a 38.6(9) keV cyclotron resonance feature in the broad 
band spectrum obtained with BeppoSAX observations from 1996. The inferred magnetic field 
strength at the neutron star surface was (3.2±0.1) xl0 12 (l+z) G, where z is the gravitational 
redshift. In addition, they assumed that the quasi-periodic oscillation frequency was due to 
the beating between the the pulse frequency and the Keplerian motion at the magnetospheric 
radius. Assuming a distance 5<dfc pc <13, they obtained 2.4 <Bi 2 <6.3. Our measurements 
of the energy of the cyclotron line feature with Suzaku 36.8+gg and 37.7(7) are consistent 
within errors. They are in agreement with the BeppoSAX results, taking into account 
the uncertainties of the measurements. The fact that the value of the center of energy of 
the cyclotron line feature did not exhibit a notable change at different luminosity states , 


resembles the only case of the Be/X-ray binary sys 
but contrary to sources like Her X-l or V 0332+53 (Pottschmidt et aj 
therein) . 


em A 0535+26 (Caballero et al.l 120091) . 
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and references 


More recently, IZhang and Li (120101) proposed a new explanation for the 2008 torque 
reversal. They claimed that winds or outflows from the neutron star and the accretion disk 
may play an important role in accounting for the spin-down in the disk-fed neutron stars 
like 4U 1626-67. The essential idea is that the spin-down is induced by stellar and disk winds 
(outflows) that take away the angular momentum of the neutron star. Thus a significant 
decrease of the mass transfer rate and possible propeller effect are not required. The main 
limitation of this model is the mechanism originating these winds. 


In general, torque switching models explain this behavior by the effect of the accretion 
flow, and the neutron star is assumed to be acting as a solid body. Observations of glitches 
(and glitch-like events in accreting pulsars) as well as unusually rapid spin-down phases 
suggest that the this may not be the case, and the contribution of torque accumulated in 
the interior of neutron stars may be significant. 
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Fig. 1. — The main panels show the Suzaku broad-band spectrum before (left panel) and 
after (right panel) the 2008 torque reversal. Black points denote XIS data (0.5-10 keV), PIN 
data points are in red (14-100 keV) and GSO in purple (~60-100keV). The residuals after 
fitting a simple model are shown in the middle panel, with no line emission and cyclotron 
line feature included. The bottom panel shows the residuals after fitting the line complex at 
~lkeV, the ~6.5keV iron line and the ~37keV cyclotron line. 


Table 1: BEST BROAD BAND FIT SPECTRAL PARAMETERS’ 71 ^ 


Observation 

JVg 

a* & 

E* 

111 cut 

E Fold 


^BBody & 

Flux 6 

constant 

X?(DOF) 



norm 



J 

a CRF 

norm® 


factor^ 


400015010 

1.15(2) 

0.847(9) 

19.3(5) 

13(1) 

36.8^ 

0 231 +0 ' 01 
U.ZOi _ 0 009 

0.63(4) 

0.71(1) 

1.4(296) 

(MJD 53797) 


0.0163(2) 



3.5(7) 

1.7(7) 


0.5(3) 


405044010 
(MJD 55444) 

0.2(1) 

0.90(3) 

0.042(2) 

21.3(6) 

14(1) 

37.7(7) 

4.2(6) 

0.520(8) 

21.5(7) 

2.5(4) 

0.56(2) 

04.(2) 

1.2(293) 


mod (PHABS (8 GAUSSIANS+GAUSSIAN+ BLACKBODY +POW)*HIGHECUT *GABS)*CONST; k x 10 21 cm- 2 ; * 
Photon Index with normalization units in photons keV - 1 cm — 2 s — 1 at 1 keV; ' keV; a Xl0 — 4 (Lumin/10 39 ergs -1 )(d/10/cpc) — 2 ; 
b xlO - 9 erg cm -2 s — 1 (0.5-100 keV); Intercalibration factors (PIN/ XIS03 and GSO/XIS 03) with XIS 03 being a spectrum 
combination from detectors XIS 0 and XIS 3 
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Fig. 2. — The X-ray flux history of 4U 1626-67 relative to the flux measured by HEAO 1, 
in the same energy band, from previous works (Chakrabarty et al. (1997): circles; Orlandini 
et al. (1998): triangle; Krauss et al. (2007): stars; Camero-Arranz et al. (2011) unfilled 
squares) and the two Suzaku observations (blue diamonds) in the 0.7-60 keV band. 


Table 2: FIT RESULTS ON LINE EMISSION 



Expected 

Line Energy 

(keV) 

Intensity* 


EW (eV) 


Observation 

(keV) 

400015010 

405044010 

400015010 

405044010 

400015010 

405044010 

0 Hea 

0.568 

0.568(fixed) 

0.57(3) 

4+ 3 

-l 

1 5 075 +0 ' 01 

10.U/0_ ooog 

7.3 

13.2 

0 Lya 

0.653 

0.659(7) 

0.661(7) 

3.2(9) 

16.609(3) 

13.5 

15.6 

Ne Hect 

0.914 

0.928(7) 

0.929(4) 

1.4(3) 

8.3690(9) 

7.1 

9.3 

Ne Lyo- 

1.021 

1.025(2) 

1.031(2) 

3.7(2) 

38.399(9) 

22.3 

45.7 

Ne He/5 

1.084 

1.084 (fixed) 

1.11(3) 

0.2 (fixed) 

2-30 

1.2 

2.9 

Ne Ly/5 

1.210 

1.210 (fixed) 

1.25(2) 

0.2 (fixed) 

3(1) 

0.7 

4.6 

He-like Mgt 

1.34 

1.33(3) 

1.36(1) 

0.4(1) 

3.2(8) 

2.6 

4.8 

Ne Recombo 

1.362 

1.48(2) 

1.51(2) 

0.9(2) 

4.3(8) 

8.3 

14.9 

Fe K-fluor. line 

6.4 

6.3(2) 

6.73(7) 

1 ^+ 0 - 6 
110 -0.5 

-j ryj +U.U7 
1 1 -0.1 

70 

24 


* photons 10 _4 cm _2 s _1 ; a the widths of all the lines were fixed to the ASCA values (see 
text); for but the iron line the widths were not fixed (values obtained: 0.8(2) eV and 0.12(2) 
eV, respectively); 1 the energy of this line is so close to the Ne recombination continuum 
and the Mg K lines, and therefore we cannot exclude a different origin for the detected line. 






